XMM-Newton Observation of Fe Ka emission from a BAL QSO: Mrk 231 

T. J. Turner^'^ & S. B. Kraemer^ 



ABSTRACT 



for the Fe Ka line detected in another BAL QSO, H1413 + 117. 

Subject headings: galaxies: active - galaxies: individual (Mrk 231) - galaxies: nuclei - 
galaxies: Seyfert 



2 I We present results from a 20 ksec XMM-Newton observation of Mrk 231. EPIC 

O I spectral data reveal strong line emission due to Fe Ka, which has rarely been detected in 

^1 this class, as BAL QSOs are very faint in the X-ray band. The line energy is consistent 

^ I with an origin in neutral Fe. The width of the line is equivalent to a velocity dispersion 

I ~ 18,000 km s~^ and thus the line may be attributed to transmission and/or reflection 

from a distribution of emitting clouds. If, instead, the line originates in the accretion 
disk then the line strength and fiat X-ray continuum support some contribution from a 
^ I reflected component, although the data disfavor a model where the hard X-ray band is 

O ' purely reflected X-rays from a disk. The line parameters are similar to those obtained 
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^1 1. Introduction - the BAL QSO phenomenon 

> 

^ , The ejection of matter at moderate to high velocities appears to be a common phenomenon in Active 
^ [ Galactic Nuclei (AGN). Details of the accretion and ejection processes provide fundamental insight 
into the fueling of the AGN, and the effect of an active nucleus on its environment. Moderate 
velocity (hundreds of km s~^) outflows have been observed in many Seyfert 1 galaxies (e.g. Kaspi 
et al. 2002). Broad Absorption Line Quasars (BAL QSOs) are the high- velocity systems, showing 
blueshifted optical and UV absorption lines from resonant transitions of ionized species such as Civ, 
Siv, Nv and Ovi. Fabian (1999) suggested that BAL QSOs are in an early evolutionary phase when 
the mass of gas around a black hole is beginning to get blown away to reveal a QSO. 

The BAL absorption lines are thought to arise in radiatively-driven material along the line-of-sight, 
with outflow velocities in the range 5000-30,000 km s~^ (e.g. Turnshek et al. 1988). BAL QSOs 
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comprise ~ 8% of the high- luminosity QSOs (see Weyman 1997 for an overview). Models for BAL 
QSOs seek to explain how gas gets accelerated to the observed velocities. Suggestions have included 
a scenario where gas and radiation pressure lift material from the photosphere of the accretion disk 
and this gas is illuminated by radiation from the inner disk/corona (Murray et al. 1995). That 
nuclear photon flux then accelerates the wind outwards. As the inner edge of the wind is thought 
to be quite close to the central source, some "shielding" gas has been invoked to prevent the wind 
material becoming stripped, and difficult to drive. This gas is thought to be highly-ionized and have 
column density in the range = 10^^ — lO^'^cm"^. Proga ct al. (2000) found that a layer of shielding 
gas arises naturally in their models for BAL winds. An alternative suggestion is that the outflow 
consists of dense clouds that are magnetically confined and that have small filling factor (Arav, Li 
& Begelman 1994). 

Constraints to date on the BAL gas location and mass-loss rates have been based on the UV and 
optical data. X-rays may be the key to understanding BAL QSOs as they offer the opportunity 
to probe closer to the nucleus than studies in the optical and UV bands. Interesting X-ray results 
have recently emerged with detection of some broad absorption features from the Icnscd QSO, APM 
082769+5255 (Chartas et al. 2002), indicating an origin in gas components outflowing at ~ 0.2 — 0.4c. 
However, BAL QSOs are very faint, and detection of significant spectral features in the X-ray band 
has been rare to date. For example, while Fe Ka emission is commonly seen in Seyfcrt 1 galaxies (e.g. 
Nandra et al. 1997) and valued as an indicator of conditions very close to the black hole, this fine has 
been detected in just a few BAL QSOs, e.g. the gravitationally-lensed 'Cloverleaf QSO, H1413 -1-117 
(Oshima et al. 2001). Here we report on a brief XMM-Newton (hereafter XMM) observation of the 
BAL QSO Mrk 231, yielding a rare detection of an Fe Ka emission line. 



2. The Properties of Mrk 231 



Mrk 231 (z=0.042) is an ultraluminous infrared galaxy (ULIRG) and one of the strongest known Fell 
emitters. Furthermore, Mrk 231 has one of the highest bolometric luminosities in the local (z < 0.1) 
universe with Ls- loooMm ~ 4 X 10^^2,0 and Lbol > lO^^erg s'^ (Soifer et al. 1987). Bryant & Scoville 
(1996) find that starburst activity can only account for 40% of the IR luminosity in Mrk 231, an active 
nucleus appears to provide the rest of the energy. Another extreme property of Mrk 231 is its high 
polarization, ~ 20% at 2800 A (e.g. Smith et al. 1995), most hkely due to scattering by dust. Broad 
optical emission lines are observed similar to those defining the Seyfert 1 class. However, no narrow 
emission lines are seen except for O [ll] A3727. The absence of narrow emission lines is sometimes 
observed in quasi-stellar objects (QSOs) but rarely in Seyfert galaxies, this, and the high bolometric 
luminosity of Mrk 231 supported the idea that the nucleus of this source is similar to a QSO. This 
O [ll] emission line is consistent with production in gas which is ionized by hot young stars (Smith 
et al. 1995) rather than by ultraviolet radiation from the nucleus. This possibility seems attractive 
since the nuclear UV, extreme-UV and soft X-ray photons arc likely to be heavily attenuated by the 
dense circumnuclear absorber, i.e. the narrow-line gas may not receive ionizing radiation. 



Boroson et al. (1991), Rudy et al. (1985) and Rupke et al. (2002) are among those discussing several 
distinguishable Broad Absorption Line (BAL) systems evident in optical and UV observations of 
Mrk 231. The most prominent is defined by Nal, Can and Hel lines originating in gas with an outfiow 
velocity ~ 4200 km s"^ relative to the line emission seen in the source. Other systems appear and 
disappear between observing epochs, showing velocities from ~ 200 km s~^ infiow to ~ 8000 km s~^ 
outflow (e.g. Boroson et al. 1991; Kollatschny, Dietrich & Hagan 1992; Forster et al. 1995). Rudy et 
al. (1985) showed that the shapes of absorption features and velocities in the absorbing gas systems of 
Mrk 231 arc the same as those observed in BAL QSOs and suggested the gas is probably accelerated 
by the same mechanism. Nevertheless, the nature of the absorbing systems in Mrk 231 continued 
to be debated. The width of the CiV resonance line is used as a key indicator of the BAL nature 
of an AGN (Weyman 1997). Heavy reddening has made that fine width difficult to constrain in 
Mrk 231 until an HST spectrum settled the case in favor of a BAL QSO designation (Gallagher et al. 
2002; hereafter G02). Most BAL systems are dominated by high- ionization species, only ~ 15% have 
absorption by low-ionization systems such as MgllAA2796 and 2803, AlllAA1671, AllllAA1855 and 
1863 and CllAA1335 (Voit, Weymann & Korista 1993). In addition to these low-ionization systems, 
Mrk 231 shows even rarer absorption lines, such as NalD, HelA3889 and MglA2853. While Mrk 231 
docs show some of the high-ionization systems (G02) it is dominated by the low-ionization systems; 
the Mgll absorption in particular, places Mrk 231 into the so-called LoBAL subclass of low-ionization 
BAL QSOs. 



2.1. Previous X-ray Observations of Mrk 231 



The ROSAT High Resolution Instrument (HRI) showed extended soft X-ray emission consistent 
with the extent of the host galaxy (~ 10"; Turner 1999), confirmed using Chandra (G02). An early 
ASCA observation (Iwasawa 1999; Turner 1999) showed a fiat spectrum but with insufficient counts 
to distinguish between absorption or refiection models. In this paper we use the term "refiection" 
to refer to reprocessing of photons via Gompton scattering and fiuorescence by material which is 
optically-thick to electron scattering. We use the term "scattering" to refer to the case where the 
optical depth to electrons is < r < 1. A longer ASCA observation (Maloney & Reynolds 2000; 
hereafter MR) suggested the X-ray spectrum to be dominated by reflected plus scattered nuclear 
emission with a large absorbing column in the linc-of-sight. Thus MR suggested a Seyfert 2 type 
geometry and orientation for Mrk 231, with the accretion-disk system surrounded by a shell of 
absorbing gas Nh ~ 3 x lO^^cm"^ which covers direct and reflected/scattered components. The 
shell may block ionizing radiation from the optical narrow-line region, while allowing hard X-rays to 
penetrate. G02 find significant variability in the hard X-ray fiux, and suggested nuclear X-rays are 
hidden by a Gompton-thick cloud (Fig 12 of G02) existing on size-scales < lO^^cm, while scattered 
X-rays from multiple lines-of-sight reach the observer. In the context of the MR model the hard 
X-rays could vary if the reprocessor were suitably small. This may be reasonable for a source like 
Mrk 231, radiating at a large fraction of the Eddington luminosity, as the inner disk may puff up 
into a funnel (Rees 1984) yielding a reflector surface consistent with the size implied by variability 
in this case. Smith et al. (1995) suggest there may be a superposition of two polarized components 



in Mrk 231, supporting the importance of scattering in this system. 



3. 



The XMM-Newton Observation 



3.1. 



Data Reduction 



An XMM observation of Mrk 231 was performed starting 2001 June 7 UT 13:10:24 with an exposure 
time of ~ 22 ks. Here we present analysis of the archived XMM data from that observation. EPIC 
data were processed using SASS 5.4. EPIC pn and MOS data utihzed the medium filter and the 
prime full window. The EPIC data were screened to remove hot and bad pixels, and some brief 
periods of high background counts. Instrument patterns 0-12 (MOS) and 0-4 (pn) were selected. 
The screening criteria resulted in an effective exposure of 21.4 ks for the MOS and pn data. The 
background-subtracted count rates were 0.0391±0.0014 (MOSl), 0.0389±0.0014 (M0S2) and 0.130± 
0.0028 (pn) over the 0.3-10.0 kcV band. Spectra and time series were extracted from a cell of radius 
~ 30" centered on the source. Background spectra were extracted from nearby, source-free regions of 
MOS and pn, the background was ~ 6% of the source count rate. Unfortunately this brief exposure 
accumulated only 400 source counts in the combined arms of the RGS data, insufficient for spectral 
analysis. 



XMM data showed Mrk 231 to have an observed flux F2-1Q kev ~ 6.2 x 10~^^erg cm~^ s~^, during 
the June 2001 observation, corresponding to an observed luminosity -L2-10 keV ~ 2 x lO^^erg s"^ 
(assuming Hq = 75 km s"^ Mpc"-^, = 0.5 throughout). The 2-8 keV flux (5 x 10~-^^erg cm~^ s""-^) 
is consistent with that reported by G02. We found no signiflcant variations in X-ray flux within this 
short XMM exposure. 

The overall spectral shape is complex. Several previous papers have discussed the difficulty in dis- 
tinguishing between models where the active nucleus is covered by a patchy and/or ionized absorber 
(Iwasawa 1999; Turner 1999) and those where the spectrum is dominated by a reflected component 
(MR). This XMM exposure is brief and the data do not allow us to determine the precise contribu- 
tions of reflected and absorbed continua to the overall spectrum. However, progress is possible as 
the pn spectrum shows two prominent emission lines (Fig 1) plus emission consistent with a blend of 
lines from the (previously suspected) thermal component peaking between 0.5 - 1 kcV (Iwasawa 1999; 
Turner 1999; Gallagher et al. 2002, MR). EPIC data show some evidence for absorption-like features 
at ~ 1.1, 1.6 and ~ 2.4 keV. The latter appears only in pn data and is probably an instrumental 
effect. For the others, the low S/N combined with the modest energy resolution leads to ambiguity 
between emission and absorption features. These may be absorption features due to Fe, but these 
data do not allow an unambiguous determination. 



3.2. 



EPIC Results 
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Fig. 1. — Left: The EPIC pn spectrum of Mrk 231 from a 20 ksec exposure with XMM. A large flux 
of lines is also evident between 0.5 - 1 keV. A sharp dip just above 2 kev does not show up in MOS 
data and is likely to be an artifact of the instrument. The dashed vertical line shows where we would 
expect the Cvi RRC feature, at the redshift of the host galaxy. The right panel shows a close-up of 
the Fe Ka line. Dashed lines denote the expected placement of a line emitted at 6.4 keV and an edge 
at 7.1 keV, observed at the redshift of the host galaxy. 

One of the most obvious features in the pn data is immediately identified with emission from the 
K-shell of Fe. The pn detects 30 counts in the hne during this exposure. The MOS spectra have far 
fewer counts and the hne is not significantly detected in those (although the data are consistent with 
the pn results). Thus we used only the pn data to derive the line parameters given below. Detection 
of an Fe line in the EPIC spectrum from Mrk 231 is also reported within the analysis of a sample of 
Ultra-Luminous Infrared Galaxies performed by Franceschini et al. (2003). 

First we fit the Fe Ka line relative to a 'local' continuum parameterized as a powerlaw between 
4-9 keV, absorbed by the galactic line-of-sight column [Nh = 1.27 x 10^° cm~^). This yielded 
r ~ 1.0llQ 3g, indicative of complex absorption and/or a contribution from reflection, as previously 
found for this source. We measured an observed-energy for the emission line E = 6.23 ± 0.12 keV, 
corresponding to a rest-energy E = 6.49 ± 0.12 keV correcting for the systemic velocity of the host 
galaxy. The line width was a = O.IT^q'H keV with flux n = 3.85lfgg x 10^^ photons cm~^ s~^ 
and the luminosity in the line was LpeK = 1.3 x 10^^ erg s^^. The lower limit of a = 60 eV means 
the Fe Ka line has a significant breadth. (Errors are 90% confidence.) The line equivalent width is 
~ 450I204 relative to the local continuum. For comparison the spectral resolution of the pn at 6.4 
keV is cr ~ 50 eV. Fig 1 shows the data compared to lines denoting the expected observed energies 
of an emission line from neutral Fe, along with the energy of the Fe K edge which would arise in the 
same gas, both assumed to be observed at the redshift of the host galaxy. While a dip is evident at 
the energy where an edge is expected, addition of an edge did not improve the fit significantly. 



For comparison, we fit the line relative to a "pure reflection" model for data in the 2-10 keV band 
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Fig. 2. — The left-hand panel the Fe Ka line in the combined instruments provided by ASCA data, 
the dashed vertical lines represent the energies of a line at 6.4 keV and an edge at 7.1 keV observed at 
the redshift of the host galaxy. The right-hand panel shows data in the same regime from our analysis 
of the archived ACIS observation of 2000 October. The Fe Ka line is evident in both datasets, and 
the flux is consistent with that found by the pn. 

(as suggested by MR). This fit utilized the pexrav model (Magdziarz & Zdziarski 1995) from XSPEC. 
Fits were attempted assuming disk inclination angles of 45° and 0° (face-on). The (unseen) powerlaw 
was assumed to be F = 2.0 and to turn-over at 100 keV. Fits using the former disk orientation gave 
= 103/97 dof and a line equivalent width EW = 202^1^1 eV relative to the reflected continuum, 
while a face-on disk yielded = 103/97 dof, EW = 139^^3g. The best-fit equivalent width in either 
case is lower than that found relative to the absorbed powerlaw because the edge in the reflected 
continuum contributes to the curvature just above 6.4 keV, and so fewer photons are attributed to the 
line. In these two (pure reflection) disk models the best-fit line strength is about 10% of the strength 
expected, relative to the reflected continuum (George & Fabian 1991). Thus the line strength is low 
compared to the model in this picture. As the discrepancy is in the ratio of reflected line to reflected 
continuum, we can extend this result to a range of geometries. The line is too weak relative to the 
reflected continuum if reflection occurs from either a disk or a distribution of clouds. 

Alternative fits were attempted, which utilized both transmitted and reflected components, with an 
unconstrained absorber covering both. First we considered the face-on disk. We fixed the reflection 
to that expected from a flat disk illuminated by an isotropic source (i? = 1 in pexrav). This yielded 
= 93/97 (lof with a flat photon index F = 1.14+°;^^, Nh = 2.69^?;^^ x 10^^ cm'^, indicating 
additional, but undeterminable complexity to the absorber. We obtained similar results for a disk 
at 45°, with = 93/97 dof, F = 1.14+°;^^, A^^ = 2.71^?:^^ x 10^^ cm'^. Fixing, instead, the photon 
index at F = 2 yielded i? ~ 10 for both disk orientations, equivalent to the 'pure reflection' solutions, 
above. 



To confirm the presence of the Fe line, and look for flux variability we reduced the long archival 



observation by ASCA, from 1999 November 10-12. Data reduction followed the methods outlined 
for the Tartarus database (Turner et al. 2001). The ASCA spectral data were fit using all four 
instruments, and the spectrum confirms the presence of the Fe line (Fig 2). The spectrum was 
modeled constraining photon index, absorption, line energy and line width to match the values 
derived from the XMM data. This produced a good fit, with continuum flux F2-1Q kev ~ 7.5 x 
10"^^erg cm"^ s"^ and hue flux n = 5.76^2^1 x 10"^ photons cm"^ s"^ 

We also reduced the archived ACIS data from 2000 October 19-20, originally reported by G02 (Fig 2). 
Extracting ACIS CCD spectra of Mrk 231 we fit the data using the same continuum and line model 
used for the XMM and ASCA spectra, leaving only line and continuum normalizations free. This 
model yielded a line flux n — 2.48^2^1 ^ 10"^ photons cm~^ s~^. The differences in hue flux are not 
significant and so the data show the line flux to be consistent with a constant value over a baseline 
of 1.5 years. 

Turning to the soft X-ray band. Previous PSPC (Turner 1999) and Chandra ACIS data (G02) 
allowed limited spatially-resolved spectroscopy which provided two-temperature thermal models for 
the extended soft emission. While previous models account for the approximate flux in the soft 
band, they fail to produce some observed features such as the relatively isolated peak evident at 
0.471 ± 0.03 keV (Fig 1). Correcting for the redshift of the host galaxy we obtain a rest-energy 
0.490 keV, nominally consistent with Cvi RRC 0.492 keV, although most of the hue flux hes slightly 
redward of this energy (Fig 1). Conflrmation of Cvi RRC would be interesting as the width would 
constrain temperature and excitation mechanism for the gas. 



4. Discussion 

A brief XMM observation of Mrk 231 has revealed a rare detection of Fe Ka from a BAL QSO. 
The presence of strong Fe Ka emission dispels any remaining doubt as to the presence of an active 
nucleus at the heart of Mrk 231. The existence of this line is conflrmed by analysis of archived ASCA 
and ACIS data, and the line flux is consistent with a constant value between all datasets, covering 
a baseline of ~ 1.5 years. The luminosity in the line is Lpe ^ 10^^ erg s~^ which is typical of the Fe 
Ka luminosities observed in Seyfert 1 and 2 galaxies (e.g. Nandra et al. 1997; Turner et al. 1997). 
Assuming a nuclear system with an intrinsic luminosity in the Seyfert 1 range, the luminosity of 
the line would disfavor any model where we are seeing the line via scattering from optically-thin 
gas, which typically scatters only a few percent of nuclear emission into our line-of-sight in Seyfert 
2 galaxies. The other alternative is that we see only a few percent of the line flux, which happens 
to come out to a line luminosity typical of Seyfert 1 galaxies. In this case the nuclear luminosity 
would just have to be a factor 50 or so higher than a Seyfert nucleus. Different analyses have yielded 
a wide range of inferred nuclear luminosities for Mrk 231, including an estimate at 2 x lO'^'^crg s~^ 
(MR) which encompasses the latter possibility. The high equivalent width of the line, and the flat 
continuum shape are reminiscent of the properties of H1413 -|- 117 (Oshima et al. 2001). A line of 
this equivalent width could be produced from shell of gas Nh ~ 3 x lO^^cm"^ (Leahy & Creighton 



1993). This column is consistent with that found when the X-ray spectrum is fit with a powerlaw, 
partially-covered by cool gas (Turner 1999). The line is too strong to be produced by scattering from 
a shell with Nh ~ few x lO^^cm"^, as suggested to cover the central regions (MR). 

We turn now to the line energy. Wc consider the possibility that the line arises from neutral Fe. If so, 
there is no systematic velocity shift for the emitting gas. While the data disfavor models where we 
preferentially observe one side of a flow system, the line breadth is equivalent to a velocity-width of 
18500^12000 The geometry of the emitting gas remains unconstrained, it could be anything 

from a spherical distribution or ring of clouds to the conical BAL flow, so long as we see contributions 
from both sides (as the hne peak is not shifted). A specific model which fits this general picture is 
that of Arav, Li & Begelman (1994), where the line would arise in dense clouds with a small filling 
factor. The Fe Ka emission in H1413 + 117 (Oshima et al. 2001; Gallagher et al. 2002b) was also 
found to be consistent with an origin in neutral gas. In that case, emission from the BAL fiow was 
used to explain the line, with the width suggesting a fiow velocity of 12,000 km s~^. 

An alternative is that the line is produced in the accretion disk and that the broadening can be 
attributed to strong gravity, as it is in the spectra of Seyfert type AGN, where there are no high- 
velocity BAL fiows. Of course the hne width observed here is rather modest compared to the 
distortion observed in Seyfert profiles. The equivalent width of the hne is ~ 450 eV, about twice 
that expected from a source isotropically illuminating a fiat disk (George & Fabian 1991), assuming 
we see both the primary and the refiected components. This may infer the reprocessor sees more of 
the illuminating flux than the observer does, or that there is a lag between primary and reprocessed 
spectra. An alternative explanation of the line strength is that the Fe abundance is very high in 
BAL QSOs, as expected if these are young systems (Hamann & Ferland 1999). 

However, fits assuming we see only the refiected spectrum in the 2-10 keV band (i.e. the primary 
continuum is hidden) show the Fe line has just 10% of the strength expected, relative to the refiected 
continuum. This supports the case for some contribution from the primary continuum in the hard 
X-ray regime and thus disfavors pure refiection models, whether refiection occurs from a disk or from 
a distribution of clouds. 

Of course, the breadth of the line could be due simply to Compton scattering of line photons in 
material along the line-of-sight, as discussed during the early days of Fe Ka observations from 
Seyfert 1 galaxies (e.g. Mushotzky et al. 1995). If so, the lack of shift to the peak indicates the 
opacity of that gas is low, while the upper limit on line width, a, constrains the gas to kT < 0.5 keV 
(using AE/E = y/2kT/m^), following Mushotzky et al. 1995). 

Finally, it is impossible to rule out emission from a higher energy, redshifted down to 6.49 keV. A 
velocity of 24,000 km s^^ is possible if the line is from recombination in H-likc Fe. However, if the 
gas were highly ionized wc would most likely sec emission from both Hc-likc and H-likc Fe (and the 
observed feature would have to be a blend of these two). Further, it seems unlikely H1413 + 117 and 
Mrk 231 would both show single broad lines consistent with 6.4 keV by this conspiracy. Thus we 
favor an origin in neutral gas. 
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